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Abstract 

This paper addresses the need to understand the physics and chemistry involved in propellant combustion 
processes in micro-scale combustors for propulsion systems on micro-spacecraft. These spacecraft are 
planned to have a mass less than 50 kilograms with attitude control estimated to be in the 10 milli- 
Newton thrust class. These combustors are anticipated to be manufactured using Micro Electrical 
Mechanical Systems (MEMS) technology and are expected to have diameters approaching the quenching 
diameter of the propellants. Combustors of this size are expected to benefit significantly from surface 
catalysis processes. Miniature flame tube apparatus is chosen for this study because microtubes can be 
easily fabricated from known catalyst materials and their simplicity in geometry can be used in 
fundamental simulations for validation purposes. Experimentally, we investigated the role of catalytically 
active surfaces within 0.4 and 0.8 mm internal diameter microtubes, with special emphases on ignition 
processes in fuel rich gaseous hydrogen and gaseous oxygen. Flame thickness and reaction zone thickness 
calculations predict that the diameters of our test apparatus are below the quenching diameter of the 
propellants in sub-atmospheric tests. Temperature and pressure rise in resistively heated platinum and 
palladium microtubes was used as an indication of exothermic reactions. Specific data on mass flow 
versus preheat temperature required to achieve ignition are presented. With a plug flow model, the 
experimental conditions were simulated with detailed gas-phase chemistry, thermodynamic properties, 
and surface kinetics. Computational results generally support the experimental findings, but suggest an 
experimental mapping of the exit temperature and composition is needed. 

Introduction 

The need for micro-scale propulsion devices is 
apparent due to the growing desire to reduce the 
size and cost of future spacecraft. These spacecraft 
are planned to have a mass of less than 
50 kilograms and require an attitude control 
estimated to be in the 10 milli-Newton thrust class 
[1]. NASA Goddard Space Flight Center is 
currently conducting feasibility studies on 
formation flying of such spacecraft and the U.S. 

Air Force is planning on developing 10 to 20 kg 
class of spacecraft, defined as “nano-spacecreaft,” 
for use in satellite constellations [2]. Chemical 
propulsion may have a distinct advantage because 
hydrocarbon energy densities are -100 times 
greater than current batteries [3] and the fact that a 


small amount of electrical energy will be budgeted 
for propulsion purposes. In order to develop 
chemical based propulsion for this class of small 
spacecraft, micro-combustion must be realized. 
These micro-combustion propulsion devices can be 
used both for primary thrust and orbit or attitude 
control. When grouped into arrays for larger thrust 
applications, they provide for unprecedented 
propulsive flexibility with a range of thrust levels 
and steady low-g acceleration. The redundancy of 
the thrusters and the potential for a spark-less 
ignition can also provide unmatched reliability. In 
a more general application, they can be used as 
igniters for larger rockets. Micro-combustors for 
power generation currently being researched are 


* Aerospace Engineer, On-Board Propulsion Branch, StevenJ.Schneider@grc.nasa.gov , Associate Fellow AIAA 
t Associate Professor, Department of Mechanical and Aerospace Engineering, cjsl5@po.cwm.edu . Senior Member AIAA. 
§ Graduate Student, Department of Mechanical and Aerospace Engineering, gabl6@po.cwru.edu . Student Member AIAA. 


NAS A/TM— 2003 -2 12613 


1 



predicted to possess energy densities of 10 times 
that of an existing battery [4]. 

If combustion is employed as a viable form of 
propulsion, the devices are anticipated to be 
manufactured using Micro Electrical Mechanical 
Systems (MEMS) technology and are expected to 
have diameters approaching that of the quenching 
diameter of the propellants. In order to overcome 
the heat loss to the wall of the combustor, the 
thruster will need to take advantage of the added 
benefits of catalytic combustion on the surface of 
the device. More research on the catalytic effects 
of different materials on a wide array of common 
fuels is still needed, especially those which are a 
high priority of in-situ resource utilization (ISRU) 
technologies. The overall objective of this study is 
to explore flammability limits, combustion 
efficiency, and modeling of propellant combustion 
over heated catalysts whose geometry approaches 
or is smaller than the order of flame thickness of 
the propellant. 

Miniature, cylindrical flame tubes on the order 
of sub-millimeter in size are favorable for 
investigation in this study because of the simple 
geometry and relative ease of fabrication. In each 
case, the entire channel can be fabricated from the 
catalyst material with minimum effort. 

In the following, we will sequentially present 
the characterization of gaseous H 2 /0 2 combustion, 
an overview of catalytic combustion, experimental 
setup, and computational methodology, which is 
followed by results and discussion. 

Flame Characteristics and Quenching Diameter 

Studies of premixed combustion in a tube have 
shown that if the tube diameter is progressively 
decreased, the flame eventually cannot propagate 
[5,6]. Thus, the so-called quenching diameter is 
the critical diameter of the tube below which flame 
propagation is not possible. This is because 
diffusion of species and heat transfer from the 
flame front to the tube wall quenches the reaction. 
It is also noted that wall quenching affects not only 
flammability limits, but also ignition phenomena 
[5,6]. 

Theoretical analysis has shown that the 
quenching diameter is on the same order as 
the flame thickness [5]. The flame thickness 
represents the overall non-equilibrium processes 
of reaction and diffusion in the flame, and hence 
provides an indication of the residence time within 
the flame. Since the flame thickness is inversely 


proportional to the burning rate, it decreases with 
increasing burning intensity. As such, near- 
stoichiometric flames would have a relatively 
smaller flame thickness as well as quenching 
diameter due to stronger burning intensity. 

Spatially resolved profiles of temperature and 
heat release rate for a freely-propagating H 2 /0 2 
flame in the flame coordinate are computed using 
PREMIX code [7] with detailed chemistry, 
thermodynamic properties, and transport 
properties. The H 2 /0 2 kinetic mechanism is taken 
from Mueller et al. [8]. Calculations were run for 
freely propagating H 2 /0 2 premixed flames under 
adiabatic and non-adiabatic conditions. 

Flame dimensions can be characterized using 
various methods of analysis. In the first method, a 
flame thickness, 8, can be obtained by calculating 
the difference between the final and initial 
temperatures and dividing by the maximum value 
of the temperature gradient. A plot of H 2 /0 2 flame 
thickness as a function of equivalence ratio for 
atmospheric pressure and pressure of 34.5 kPa 
(approximately 5 psia) is given in fig. 1. As 
expected, 8 increases with increasing equivalence 
ratio. It is also seen that flame thickness increases 
with decreasing pressure and the calculated flame 
thickness at the test pressure of 34.5 kPa and an 
equivalence ratio of 10 (mixture ratio MR = 0.8) is 
4 mm, which is greater than the diameter of either 
the 0.4 mm platinum or the 0.8 mm palladium 
microtubes used in this catalytic micro-combustion 
experiment. 



Figure 1 . Plot of flame thickness and reaction zone thickness 
versus equivalence ratio for p = 1 atm and 5 psia. 

In the second method, we define a reaction 
zone as the full width of the half maximum of 
the heat release profile. Figure 1 also shows the 
corresponding reaction zone thickness. At the 
equivalence ratio of 10 (MR = 0.8), the reaction 
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zone thickness is ~1 mm at atmospheric pressure 
and initial temperature of 300 K and increases to 
~3 mm when the pressure is reduced to 34.5 kPa. 
Experimental tests in the microtubes were run 
mostly between 1.2 and 0.4 mixture ratios 
(equivalence ratios of 6.7 to 20, respectively). At 
these mixture ratios and the low pressure of the 
tests the reaction zone thicknesses are greater than 
the diameter of either the 0.4 mm platinum or the 
0.8 mm palladium microtubes used in this catalytic 
micro-combustion experiment. 

In order to determine the flammability limits 
and capture the extinction turning point, the 
PREMIX code was modified to include radiative 
heat loss and the capability of solving the 
equations at the singular points. When continuing 
the calculation from the upper (stable) branch to 
the lower (unstable) branch by passing through the 
turning point, numerical difficulty, especially with 
detailed chemistry, is encountered because the 
Jacobian matrix used in the Newton’s method 
becomes singular at the turning point. To solve this 
problem, the flame controlling continuation 
method [9] is employed in this study. Radiative 
heat loss from H 2 0 is considered since it is the 
primary radiative species in hydrogen flames. The 
corresponding Planck mean absorption coefficient 
is a function of temperature and is obtained from 
Ref. [10]. The analysis was done to test the rich 
flammability limits of the H 2 /0 2 flames. Figure 2 
shows the variation in the mass burning flux with 
equivalence ratio at varying pressures. It is seen 
that the burning flux decreases with increasing 
equivalence ratio. The turning point, beyond which 
steady flame propagation is not possible, defines 
the rich flammability limit. 



Figure 2. Mass burning flux as a function of equivalence ratio 
for H 2 /0 2 flames at various pressures. The turning point of the 
response curve defines the rich flammability limit. 


Calculations also show that the rich limit of the 
H 2 /0 2 flame increases with decreasing pressure. 
We further note that for near rich limit hydrogen 
flames, in which the Lewis number is greater than 
unity, pulsating instability has been observed to 
take place in mixtures with fuel concentrations 
slightly below that of the rich limit [11-13] and the 
flame is found to oscillatorily extinguish well 
ahead of the static extinction limit [12,13]. 
Therefore, the flammability range of a H 2 /0 2 
mixture is expected to be narrowed when the 
intrinsic oscillatory nature of flame propagation is 
taken into account. 

Overview of Catalytic Combustion 

The principal advantage of using catalysts for 
primary combustion reactions is their ability to 
sustain reactions that are well beyond the 
flammability limit for homogeneous gaseous 
combustion. During the last three decades, 
catalysts have been investigated for both 
augmenting primary heat release in gas turbines 
and post-combustion treatment of pollutant 
emission in aircraft, automobile, and power 
generation applications. For generating power, the 
focus is on conditions at which catalytic reactions 
are used to promote the primary combustion 
reactions. 

A heterogeneous catalyst typically adsorbs the 
reactants and subsequently promotes critical 
reactions at certain adsorption sites. Surface 
reactions occur and then the products are desorbed 
and along with heat are transported into the flow 
by diffusion, thereby increasing the overall 
reaction rate. 

The reacting flow processes within the 
catalytic bed are complicated. Take a premixed 
fuel/air mixture flowing through a long channel 
with catalytic surfaces as an example [14]. Near 
the inlet section, the fuel conversion is controlled 
by the heterogeneous reactions at the surface of 
the catalyst. In this region, the bulk gas 
temperatures are typically too low to support 
homogeneous combustion. After a certain distance 
downstream of the inlet, the heterogeneous 
reaction rates increase with increasing 
temperatures of both near-wall gas and substrate. 
As the reaction rate increases to an extent, the fuel 
conversion is no longer controlled by the surface 
reaction rate. Instead, the fuel consumption is 
now controlled by the rate at which new gas phase 
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reactants are transported to the surface, and the 
catalytic reactions at the wall will be limited solely 
by the rate of mass transfer. Further downstream 
and beyond this diffusion-limited region, thermal 
runaway of homogeneous or gas phase reactions 
occur as heat continues to diffuse into the bulk gas 
flow. These homogeneous reactions are typically 
initiated at bulk gas temperatures in the range of 
1200-1300 K. 

Based on the above description, it is 
imperative to identify how the catalyzed processes 
scale with the combustor size [15]. In the surface 
kinetics controlled region, the process is primarily 
dependent on the catalytic Damkohler number, 
which is the ratio of system residence time to the 
characteristic heterogeneous reaction time. For a 
given residence time, the catalytic Damkohler 
number is inversely proportional to the combustor 
size due to the dependence of the heterogeneous 
reaction rate on surface area [15]. The transition 
from surface reaction- to diffusion-limited 
processes and the progression of the latter are 
expected to scale with the mass transport number 
that is the ratio of the characteristic residence time 
of the gas in the reactor to the mass/heat transport 
time. Since both the heat and mass transfer times 
scale inversely with the hydraulic diameter, both 
transport processes are favored when the 
combustor dimension is reduced, as is the overall 
progression to the homogeneous reaction region 
[15]. 

Therefore, although the increase of the surface 
to volume ratio of the combustor is problematic 
for gas phase combustion, it favors catalytic 
combustion. Despite the catalytic reaction being 
generally slower than the gas phase reaction, due 
to the relative increase of surface area and the 
lower temperature of the catalytic reaction, micro- 
combustors utilizing catalytic combustion may be 
easier to implement than those utilizing gas phase 
reactions. In fact, it has been demonstrated [16] 
using a heat-recirculating engine that with a 
catalyst located in the combustion volume 
hydrogen/air mixtures can exhibit thermal 
runaway at room temperature and butane (or 
propane) can self-start at less than 500 K. 

Many unique technology issues are related to 
MEMS -based micro-propulsion systems. They 
include overall issues such as materials selection, 
fabrication, sealing, assembly, pressurization, 
testing, and characterization, as well as the 


combustion related issues. In view of the above 
discussion the combustor remains one of the most 
challenging components to miniaturize. 
Therefore, the fundamental understanding of 
various micro-combustion phenomena constitutes 
the long range goal of the present study. 

As a first step, the present microtube 
experiments are simulated based on PLUG code 
[17] with detailed gas -phase and surface 
chemistry. In this plug-flow model, it is assumed 
that there is no mixing in the flow direction, while 
perfect mixing in the transverse directions. 
Moreover, axial diffusion of any quantity is 
negligible relative to the corresponding convective 
term, which is particularly valid in the present 
high-speed flow conditions. The computational 
specifications will be detailed in a later section. 

We further note that Raja et al [18] recently 
carried out a critical evaluation of Navier-Stokes, 
boundary-layer, and plug-flow models for reacting 
flows in a catalytic channel. The range of validity 
of various approximations was identified through 
the comparison of the solutions obtained by using 
three different formulations. While plug-flow 
models are computationally inexpensive, the range 
of validity was found to be limited [18]. We will 
show in due course that the plug-flow 
approximation is valid in most of the conditions 
investigated herein, by comparing the competing 
time scales following the study of Raja et al [18]. 

Experimental Setup 

A test stand was arranged so that the flow path 
of premixed combustible gases was directed inside 
a microtube, fabricated from catalytic material. 
H 2 /0 2 mixtures were used in conjunction with a 

0. 4 mm I.D. platinum microtube and a 0.8 mm 

1. D. palladium microtube. Miniature flame tubes 
were chosen for experimentation because of their 
simplistic geometry that can be used in 
fundamental simulations as well relative ease of 
fabrication. The tubes will also be resistively 
heated with different values of heat flux using a 
DC power supply. The 0.4 mm I.D. platinum tube 
has a length of 100 mm, of which 70 mm was 
resistively heated, and a wall thickness of 
0.15 mm. The 0.8 mm I.D. palladium tube has a 
190 mm length, of which 160 mm was resistively 
heated, and a 0.1 mm wall thickness. In both cases, 
the middle of the tube was resistively heated, 
leaving the ends on either side to be unheated, with 
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15 mm on each side. Thermocouples were spot 
welded at locations of 20, 50, and 80% of full 
length on the outer surface of the heated section of 
the tubes. The entire test apparatus consists of the 
microtube, a mass flow control system, electrical 
power supply, and data acquisition system. 
Schematics and photograph are shown in figs. 3 
and 4 respectively. The test apparatus is then 
located inside a vacuum test cell located at NASA 
Glenn Research Center, which can maintain an 
ambient pressure of 1.38 kPa (10 torr). 


Computational Methodology 

Figure 5 shows the microtube configuration 
used in the simulations, where U is the inlet 
velocity, L the tube length, and D the tube 
diameter. The plug flow reactor is modeled to be 
an empty tube with a negligible wall thickness. 
Part of the tube is resistively heated with uniform 
heat flux Q, while the first section of length l x and 
the end section of length l 2 are not heated. 
Governing equations for solving the steady-state 
plug flow reactor is listed below. 
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Figure 3. Schematic of microtube test articles. 



Figure 4. Microtube test apparatus (0.8 mm I.D.). 


Low mixture ratios (high equivalence ratios) 
were tested to explore the fuel rich flammability 
limits of the propellants. Fuel rich conditions were 
selected to avoid melting of the catalyst and 
oxidation loss of catalyst material. Each of these 
chemicals is non-toxic and can be produced in 
orbit by water electrolysis [19]. Detailed gas 
chemistry, thermodynamic properties, transport 
properties, and surface kinetics for H 2 /0 2 /Pt and 
H 2 /0 2 /Pd systems are also available for numerical 
simulations. 
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Figure 5. Computational configuration of microtube. 


Continuity: 

. du .dp ^ . TI7 

pA— + uA — = a 2_gkW k 
ax ax gas 

Species: 

dY Kg 

puA—L+Y^g.W, = W k (g k a + d) k A) 

ax g as 

Energy: 

Jx dY k — dT du 
puA(2^ h k — + C p + u — ) + 


dx 


dx dx 


=a{Q + Y j g k W k h k ) 

gas ^ gas gas 


Momentum: 

.dp .du a 2 , 4^ • f\ 

A— + puA— + -pu f + ua2^g k W k = 0 

dX dX 1 gas 

Site Fractions: 

i, = 0. 


Here v is the axial coordinate, p the density of the 
gas mixture, u the gas velocity, A the cross- 
sectional area, a the surface area per unit length, 
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K g the total number of gas phase species, W k the 
molecular weight of species k , g k the molar 
production rate of species k due to all surface 
reactions, Y k the mass fraction of species k , Q) k the 
molar production rate of species k due to gas- 
phase reactions, h k the specific enthalpy of species 
k , C p the mean heat capacity per unit mass of the 
gas, T the gas temperature, p the pressure, / the 
friction factor, and s k the net production rate of 
surface species k. The friction factor is a function 
of the local Reynolds number based on the tube 
diameter, Re D . Specifically, for laminar flow 
f - 16/Re D , while the Blasius formula is used for 
turbulent flow. In addition, although detailed 
thermodynamic properties are considered, for 
simplicity the mixture viscosity is assumed to be a 
power dependence on temperature [20] when 
calculating the local Reynolds number. 

The underlying assumptions in the plug-flow 
model are that axial diffusion is negligible 
compared to axial convection and the radial 
transport times are much smaller than the axial 
convection time. As such, Raja et al [18] 
suggested practical bounds on range over which 
the plug-flow approximation is valid, namely 
D/L«Re D Sc«L/D, where Sc is the Schmidt 
number. Based on this criterion, for the present 
microtube configuration, the plug-flow model 
would be valid for most experimental conditions 
with tube diameter being less than 0.8 mm. 

The gas-phase reaction mechanism was taken 
from Mueller et al [8], involving eight gaseous 
species: H 2 , 0 2 , H 2 0, O, OH, H, H0 2 , and H 2 0 2 . 
The surface reaction mechanism for H 2 /0 2 
mixtures on platinum was taken from 
Duetschmann et al [21], including five surface 
species: Pt(s), H(s), H 2 0(s), OH(s), and O(s). 

Results and Discussion 

Experimental 

Experimental results consist of a determination 
of the minimum catalyst preheat temperature for 
initiating/supporting combustion in the two sub- 
millimeter diameter tubes. Data are obtained with 
the tubes resistively heated by the power supply in 
voltage control mode and with reactive premixed 
gases passed through the tubes. Tube temperature 
and inlet pressure are monitored for an indication 
of exothermic reactions and composition changes 
in the gases. 


Temperature data for the 0.8 mm I.D. 
palladium tube preheated to a peak temperature of 
477 K and showing an exothermic reaction is given 
in fig. 6. The total mass flow is 0.00123 g/sec of 
gaseous oxygen and gaseous hydrogen at an 
oxidizer to fuel mixture ratio of 0.8 (equivalence 
ratio 10.0). The thermocouples were located as 
shown in fig. 3. The resistivity of palladium at 
these temperatures is 17.2xl0“ 8 Q-m and the 
voltage drop across the tube is 1.21 Volts, giving 
an estimated power input of 15 Watts. At 
10 seconds into the test TCI has dropped to 338 K, 
TC2 has dropped to 412 K as the gas picks up heat, 
but TC3 has risen to 51 1 K as the gases are heated 
by catalyzed reactions between TC2 and TC3. At 
10 seconds into the test TC3 starts a rapid rise to 
750 K. Heat from this reaction is conducted 
through the tube upstream toward TC2 which starts 
a rapid temperature rise to 815 K at 16 seconds 
into the test. Further conduction to the TCI 
position is indicated by its rapid temperature rise at 
40 seconds into the test. Steady state is reached at 
TC3 at approximately 60 seconds into the test. 
The test was terminated at 100 seconds, however, 
before steady state was reached at TCI. 


0.8 mm Palladium Tube, MR= 0.8, 0.00123 g/s 



Figure 6. Temperature data for 0.8 mm I.D. palladium tube, 
MR = 0.8, 0.00123 g/sec. 

The pressure in a plenum upstream of the tube 
is given in fig. 7. The pressure at the start of the 
test is 2.7 kPa in the altitude test facility. The 
pressure rises rapidly to 34 kPa at 3 seconds into 
the test and stays there until about 10 seconds into 
the test when catalyzed reactions start to drive the 
pressure up. The pressure then rises to 59 kPa at 
71 seconds into the test and stays there until the 
test is terminated at 100 seconds. This increase in 
the inlet pressure suggests that the flow is choked 
within the catalyst tube due to a rise in viscosity 
and friction factor as the gases are heated by 
the exothermic reaction. Subsequent plug-flow 
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calculations also demonstrated that the gas flow in 
the microtube is choked when the imposed heat 
flux exceeds a critical value. 



Figure 7. Pressure data for 0.8 mm I.D. palladium tube, 
MR = 0.8, 0.00123 g/sec. 



Figure 8. Non-ignition temperature data for 0.8 mm I.D. 
palladium tube, MR = 0.8, 0.00123 g/sec. 



Time - sec 


Figure 9. Non-ignition pressure data for 0.8 mm I.D. 
palladium tube, MR = 0.8, 0.00123 g/sec. 

Temperature data for the 0.8 mm I.D. tube 
preheated to a peak temperature of 422 K and 
showing no reaction within the heated section 
are given in fig. 8. The total mass flow is 
0.00123 g/sec of gaseous oxygen and gaseous 
hydrogen at an oxidizer to fuel mixture ratio of 0.8. 
Some reactions between TC2 and TC3 cause the 
temperature of TC3 to rise to a peak of 430 K at 
6 seconds into the test, otherwise all of the 
temperatures fall to a steady state profile of TCI at 
321 K, TC2 at 351 K, and TC3 at 377 K. The 
corresponding pressure in a plenum upstream of 


the tube is given in fig. 9. The pressure rises 
rapidly and reaches a pressure of 33 kPa at 
15 seconds into the test. 

Temperature data for the 0.8 mm I.D. tube 
preheated to a peak temperature of 589 K at a 
mixture ratio of 0.6 (equivalence ratio 13.3) and a 
mass flow of 0.00160 g/s is given for a reacting 
case in fig. 10. The corresponding pressure data is 
given in fig. 11. This equivalence ratio is above the 
rich flammability limit predicted in fig. 2. If the 
preheat temperature is reduced to approximately 
530 K, the temperature data shown in fig. 12 
indicate that the reactions have not occurred within 
the heated section. 



Time ~ sec 

Figure 10. Temperature data for 0.8 mm I.D. palladium tube, 
MR = 0.6, 0.00160 g/sec. 



Time ~ sec 

Figure 11. Pressure data for 0.8 mm I.D. palladium tube, 
MR = 0.6, 0.00160 g/sec. 



TCI 

TC2 

TC3 


Time ~ sec 

Figure 12. Non-ignition temperature data for 0.8 mm I.D. 
palladium tube, MR = 0.6, 0.00160 g/sec. 

Temperature data for the 0.4 mm I.D. platinum 
tube preheated to a peak temperature of 640 K 
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which shows an exothermic reaction is given in 
fig. 13. The total mass flow is 0.000412 g/sec of 
gaseous oxygen and gaseous hydrogen at an 
oxidizer to fuel mixture ratio of 0.8 (equivalence 
ratio 10.0). The resistivity of platinum at these 
temperatures is 23.4xl0“ 8 Q-m and the voltage 
drop across the tube is 1.41 Volts, giving an 
estimated power input of 31 Watts. At 30 seconds 
into the test, the thermocouples have dropped to 
their minimum values of TCI at 400 K, TC2 at 
480 K and TC3 at 550 K. The remaining 

30 seconds of the test show an exothermic reaction 
taking place giving a rise in temperature at TC3 to 
650 K at the end of the test. This test was 
terminated before steady state was reached. The 
pressure in a plenum upstream of the tube is given 
in fig. 14. The pressure rises to 43 kPa at 
30 seconds into the test and then rises to 48 kPa at 
the end of the test, suggesting a choked flow due to 
a rise in viscosity and heat lease through the 
exothermic reaction. 
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Figure 13. Temperature data for 0.4 mm I.D. platinum tube, 
MR = 0.8, 0.000412 g/sec. 
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■TC3 



Figure 14. Pressure data for 0.4 mm I.D. platinum tube, 
MR = 0.8, 0.000412 g/sec. 

Temperature data for the 0.4 mm I.D. tube 
preheated to a peak temperature of 585 K which 
shows only a very mild reaction are given in fig. 
15. The total mass flow is 0.000412 g/sec of 
gaseous oxygen and gaseous hydrogen at an 
oxidizer to fuel mixture ratio of 0.8. Convection 


between TC2 and TC3 causes the temperature at 
TC3 to rise to a peak of 584 K at 10 seconds into 
the test, otherwise all of the temperatures fall to a 
steady state profile of TCI at 380 K, TC2 at 450 K 
and TC3 at 501 K at 25 seconds into the test and 
remain these levels except for a slight rise in TC3 
until the test is terminated at 69 seconds. The 
pressure in a plenum upstream of the tube is given 
in fig. 16 for this case. The pressure at the start of 
the test is 2.5 kPa and rises rapidly to 43 kPa at 
30 seconds into the test as it did in the reacting 
case shown in fig. 14. The pressure then rises to a 
plateau of 47 kPa at 65 seconds into the test 
indicating a mild viscosity change and remains 
there until the test is terminated at 69 seconds. 



Figure 15. Non-ignition temperature data for 0.4 mm I.D. 
platinum tube, MR = 0.8, 0.000412 g/sec. 



Figure 16. Non-ignition pressure data for 0.4 mm I.D. 
platinum tube, MR = 0.8, 0.000412 g/sec. 

Temperature data for the 0.4 mm I.D. tube at a 
mixture ratio of 0.6 (equivalence ratio 13.3) and a 
mass flow of 0.000537 g/sec is given for a reacting 
case in fig. 17. The peak preheat temperature of 
640 K was used. Note that this equivalence ratio 
is above the rich flammability limit predicted in 
fig. 2. If the peak preheat temperature is reduced to 
585 K, the temperature data shown in fig. 18 
indicate that less reactions have occurred. 
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Figure 17. Temperature data for 0.4 mm I.D. platinum tube, 
MR=0.6, 0.000537 g/sec. 



TCI 

TC2 

TC3 


Time ~ sec 

Figure 18. Temperature data for 0.4 mm I.D. platinum tube, 
MR = 0.6, 0.000537 g/sec. 



Total Mass Flow - g/sec 

Figure 19. Ignition data for 0.8 mm I.D. palladium tube. Open 
symbols are non-ignition within the heated section and solid 
symbols are ignition as indicated by tube temperature rise. 

These data and data at other mass flows and 
fuel rich mixture ratios are shown on mass flow 
versus preheat temperature plots in fig. 19 for the 
0.8 mm I.D. palladium tube and fig. 20 for the 
0.4 mm platinum tube. Using temperature and 
pressure rise as an indication of reactions the solid 
symbols indicate those test conditions where 
reactions within the heated section were observed. 
As mass flow decreases the preheat temperature 
required to initiate reactions decreases. For the 
0.8 mm palladium tube the lowest preheat 


temperature for reactions was 340 K at a mixture 
ratio of 0.8 and a mass flow of 0.000128 g/sec. 
Above a mass flow of 0.0022 g/sec reactions did 
not occur at any preheat temperature in the 0.8 mm 
I.D. palladium tube, 190 mm long with 160 mm 
resistively heated. 

For the 0.4 mm I.D. platinum tube, the lowest 
preheat temperature for reactions was 480 K 
at a mixture ratio of 1.0 and a mass flow of 
0.000187 g/sec. Above a mass flow of 
0.00054 g/sec reactions did not occur at any 
preheat temperature in the 0.4 mm I.D. platinum 
tube, 100 mm long with 70 mm resistively heated. 
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Figure 20. Ignition data for 0.4 mm I.D. platinum tube. Open 
symbols are non-ignition within the heated section and solid 
symbols are ignition as indicated by tube temperature rise. 


Computational 

While a detailed transient simulation on 
microtube experiments will be carried out in the 
future, as a first step the present steady-state plug- 
flow model would provide insights into the critical 
conditions leading to catalytic ignition as well as 
the axial variations of gas temperature, flow 
velocity, gas pressure, gas-phase composition, and 
site fraction of surface species. Platinum tube is 
first studied. The inlet mixture temperature and 
the oxidizer to fuel mixture ratio are kept constant 
at 300 K and 1.0, respectively. For this rich H 2 /0 2 
mixture, the inlet viscosity is 0.12 millipoise. 
Moreover, as demonstrated in fig. 1, the diameters 
of the platinum microtubes studied herein are 
smaller than the characteristic flame thickness. In 
the following, effects of imposed heat flux (Q), 
inlet velocity (U), and micro tube diameter on 
catalytic ignition are systematically examined. 


Effect of Input Heat Flux 

Consider a platinum tube of D=0.4 mm, 
L = 100 mm, and = 15 mm. With inlet 
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pressure of 0.5 atm and inlet velocity of 50 m/s, the 
inlet Reynolds number is 128. Figure 21 plots the 
axial variations of gas temperature and pressure 
from the inlet with various input heat fluxes. The 
microtube is heated within the range bounded by 
two long dashed lines. It is seen that with Q = 0.5 
W/cm 2 , catalytic ignition cannot be achieved for 
the given tube length. With increasing Q to 0.55 
W/cm 2 , ignition occurs near the exit of microtube 
and outside the heated section. Note that 
experimentally tube temperatures were monitored 
for an indication of catalytic ignition. Since 
thermocouples were spot welded within the heated 
section, computed results of fig. 21 imply that the 
experimentally-determined critical heat flux 
leading to ignition may be overestimated. 
Therefore, measurements of exit temperature and 
composition are needed for accurate determination 
of the ignition condition. 


Inlet Conditions: 300 K, 0.5 atm, 50 m/s 



Figure 21. Axial profiles of temperature and pressure in a 
0.4 mm I.D platinum tube with various input heat fluxes. Inlet 
conditions: MR = 1.0, 300 K, p = 0.5 atm, and U = 50 m/s. 
The section bounded by two long-dashed lines is heated with 
uniform heat flux. 

Figure 22 demonstrates the variations of Mach 
number corresponding to fig. 21. For a given 
heat input, Mach number increases towards the 
tube exit, while pressure decreases as the flow 
accelerates. Figure 22 also shows that further 
increasing the imposed heat flux beyond 
0.57515 W/cm 2 , the flow is choked somewhere in 
the tube. Once the flow is choked, the steady state 
solutions fail to exist for the given inlet conditions. 

For Q = 0.57515 W/cm 2 , fig. 23 further plots 
the spatially-resolved profiles of the mole fractions 
of gas species, including H 2 , 0 2 , and H 2 0, as well 
as the surface site fractions of Pt(s) and H(s). It is 
seen that for the present fuel rich condition the 
surface is almost covered entirely by H(s) prior to 


the initiation of significant exothermic reactions. 
With increasing temperature, the surface coverage 
of Pt(s) increases and eventually dominates, 
leading to the consumptions of H 2 and 0 2 and the 
production of H 2 0. 


Inlet Conditions: 300 K, 0.5 atm, 50 m/s 



Figure 22. Axial profiles of temperature and Mach number in 
a 0.4 mm I.D platinum tube with various input heat fluxes. 
Inlet conditions: MR=1.0, 300 K, p=0.5 atm, and U=50 m/s. 
The section bounded by two long-dashed lines is heated with 
uniform heat flux. 


Inlet Conditions: 300 K, 0.5 atm, 50 m/s 



Figure 23. Axial profiles of gaseous species mole fractions 
and surface coverage in a 0.4 mm I.D platinum tube with an 
input heat flux of 0.57515 W/cm 2 . Inlet conditions: MR = 1.0, 
300 K, p = 0.5 atm, and U = 50 m/s. The section bounded by 
two long-dashed lines is heated with uniform heat flux. 


Effect of Inlet Velocity 

For a given heat input of 0.15 W/cm 2 in a 
platinum tube of D = 0.4 mm, L = 100 mm, and 
t\ - I 2 - 15 mm, figs. 24 and 25 compare the axial 
distributions of temperature, pressure, and Mach 
number at three different inlet velocities, namely 
U = 10, 20, and 25 m/s. The corresponding inlet 
Reynolds numbers are 17.4, 34.8, and 43.5, 
respectively. The inlet pressure is 0.34 atm. As 
expected, lower inlet velocity facilitates ignition 
due to longer residence time. It is also evident that 
the flow is accelerated both frictionally and 
thermally (if with significant heat release). 
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Inlet Conditions: 300 K, 0.34 atm 



Figure 24. Comparison of axial profiles of temperature and 
pressure in a 0.4 mm I.D platinum tube with various inlet 
velocities. The input heat flux is 0.15 W/cm 2 . Inlet conditions: 
MR = 1.0, 300 K, p = 0.34 atm. 


Inlet Conditions: 300 K, 0.34 atm 



Figure 25. Comparison of axial profiles of Mach number in a 
0.4 mm I.D platinum tube with various inlet velocities. The 
input heat flux is 0.15 W/cm 2 . Inlet conditions: MR = 1.0, 300 
K, p = 0.34 atm. 


Inlet Conditions: 300 K, 0.34 atm 



Figure 26. Comparison of temperature evolutions as a 
function of residence time in a 0.4 mm I.D platinum 
tube with various inlet velocities. The input heat flux is 
0.15 W/cm 2 . Inlet conditions: MR = 1.0, 300 K, and 
p = 0.34 atm. 


Figure 26 further re-plots the temperature 
distributions of fig. 24 in terms of residence time. 
It is seen that catalytic ignition is possible when 
the local temperature is heated up to exceed a 
critical value. For the present condition, this 
critical temperature is found to be around 460 K, 
which is very close to the lowest preheat 
temperature, 480 K, found in the sub-atmospheric 
experiments. 

Effect of Gas-Phase Reactions 

For most of the conditions investigated herein, 
either the microtube diameter is smaller than the 
characteristic flame thickness of the premixture or 
the mixture composition is beyond the 
flammability limit. As such, it is not expected that 
gas-phase reactions account for the major heat 
release. To demonstrate the insignificance of the 
gas-phase reactions, computed results were 
compared with and without gas-phase reactions. 
Figure 27 compares a case with a platinum tube of 
D = 0.4 mm, L = 100 mm, Ix-G- 15 mm, and 
Q = 2 W/cm 2 . Inlet pressure and velocity are 
respectively atmospheric pressure and 20 m/s. In 
fig. 27, lines are the computed results obtained 
with both gas-phase and surface reactions, while 
symbols denote the solutions considering surface 
reactions only. For clarity, the symbols do not 
represent the actual grid distribution employed in 
the calculations. The closeness of the two sets of 
numerical data, including axial profiles of 
temperature and major species, indicates that for 
microtube modeling gas -phase reactions can be 
neglected. 


Inlet Conditions: 300 K, 1 atm, 20 m/s 



Figure 27. Axial profiles of temperature and major gaseous 
species in a 0.4 mm I.D platinum tube with an input heat flux 
of 2 W/cm 2 . Inlet conditions: MR = 1.0, 300 K, p = 1 atm, and 
U = 20 m/s. The symbols do not represent the actual grid 
distribution employed in the calculation. 
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Effect of Microtube Diameter 

We have demonstrated the importance of 
residence time leading to ignition. For a given 
global residence time, it is also of interest to 
compare the critical heat flux for initiating ignition 
with various diameters of microtube. In the 
following calculations, the inlet conditions are 
MR = 1.0, 300 K, and p = 0.6 atm. The global 
residence time, representing by L/U, and the tube 
length to diameter ratio, L/D, are fixed at 1/300 
and 250, respectively. Moreover, for all cases t\ - 
€2 - £ and //L = 0.15. Three tube diameters were 
compared, including D= 0.8, 0.4, and 0.1 mm. As 
a result, the corresponding inlet Reynolds numbers 
are 368, 92, and 5.76. 


Inlet Conditions: 300 K, p = 0.6 atm, 60 m/s 



Figure 28. Comparison of axial profiles of gas temperature in 
a 0.8 mm I.D platinum tube with various heat fluxes. 
L/D = 250 and inlet conditions are MR = 1.0, 300 K, p = 0.6 
atm, and U = 60 m/s. The section bounded by two long- 
dashed lines is heated with uniform heat flux. 


Inlet Conditions: 300 K, p =0.6 atm, 60 m/s 



Figure 29. Comparison of temperature evolutions as a 
function of residence time for fig. 28. 


Figures 28-33 plots the temperature profiles in 
terms of both axial coordinate and residence time. 
For each microtube, axial distributions of gas 


temperature were compared with various input heat 
fluxes. The critical heat fluxes resulting in ignition 
for D = 0.8, 0.4, and 0.1 mm are found to be 
around 0.83, 0.38, and 0.0795 W/cm 2 , respectively. 
Thus, the smaller the tube, the less the critical heat 
flux required. In addition, the critical heat flux is 
roughly proportional to the tube diameter. This can 
be explained based on the surface area to volume 
ratio, which varies inversely with tube diameter. 
As the preceding discussion of the physical 
processes and scaling of catalytic combustion 
systems, heterogeneous catalysis is favored as a 
result of increased surface-to-volume ratio. 


Inlet Conditions: 300 K, p = 0.6 atm, 30 m/s 



Figure 30. Comparison of axial profiles of gas temperature in 
a 0.4 mm I.D platinum tube with various heat fluxes. 
L/D = 250 and inlet conditions: MR = 1.0, 300 K, p = 0.6 atm, 
and U = 30 m/s. The section bounded by two long-dashed 
lines is heated with uniform heat flux. 


Inlet Conditions: 300 K, p =0.6 atm, 30 m/s 



Figure 31. Comparison of temperature evolutions as a 
function of residence time for fig. 30. 
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Met Conditions: 300 K, p = 0.6 atm, 7.5 m/s 



Distance from Met (cm) 

Figure 32. Comparison of axial profiles of gas temperature in 
a 0.1 mm I.D platinum tube with various heat fluxes. 
L/D = 250 and inlet conditions: MR = 1.0, 300 K, p = 0.6 atm, 
and U = 7.5 m/s. The section bounded by two long-dashed 
lines is heated with uniform heat flux. 


Met Conditions: 300 K, p = 0.6 atm, 7.5 m/s 



Figure 33. Comparison of temperature evolutions as a 
function of residence time for fig. 32. 


Summary 

New micro-spacecraft will require new micro- 
propulsion technology, which will require the 
development of micro-combustors. The heat 
release rate in the combustor scales with volume, 
while heat loss rate scales with surface area. 
Consequently, heat loss eventually dominates over 
heat release when the combustor size becomes 
smaller, thereby leading to flame quenching. The 
limitations imposed on chamber length and 
diameter has an immediate impact on the degree 
of miniaturization of a micro-combustor. One 
viable combustion alternative is to take advantage 
of surface catalysis. The favorable effects of 
surface catalysis associated with the high surface- 
to-volume ratio as the combustor size is reduced 
can offset the wall quenching effects. Flame 
thickness and reaction zone thickness calculations 


predict that the diameters of our test apparatus are 
below the quenching diameter of the gaseous 
oxygen and gaseous hydrogen propellants in sub- 
atmospheric tests. 

Miniature flame tube apparatus is chosen for 
study because micro-tubes can be easily fabricated 
from known catalyst materials and their simplicity 
in geometry can be used in fundamental 
simulations to more carefully characterize the 
measured heat transfer and pressure losses for 
validation purposes. Data on the temperature and 
pressure rise as a function of time in resistively 
heated platinum and palladium microtubes were 
used as an indication of exothermic reactions in a 
series of tests varying mass flow, preheat 
temperature, and fuel rich mixture ratio. As mass 
flow decreases the preheat temperature required to 
initiate reactions is shown to decrease. 

For the 0.8 mm palladium tube, 190 mm long 
with 160 mm resistively heated, the lowest preheat 
temperature for reactions was 340 K at a mixture 
ratio of 0.8 and a mass flow of 0.000128 g/sec. 
Above a mass flow of 0.0022 g/sec reactions did 
not occur at any preheat temperature in the 
0.8 mm I.D. palladium tube. For the 0.4 mm I.D. 
platinum tube, 100 mm long with 70 mm 
resistively heated, the lowest preheat temperature 
for reactions was 480 K at a mixture ratio of 1.0 
and a mass flow of 0.000187 g/sec. Above a mass 
flow of 0.00054 g/sec reactions did not occur at 
any preheat temperature in the 0.4 mm I.D. 
platinum tube. 

A plug flow model with detailed gas-phase 
and surface chemistry was used to provide 
insights into the present microtube experiments. 
Computed results for the H 2 /0 2 mixtures reacting 
on platinum show that there exists a critical heat 
flux input for critical ignition. Further increasing 
the imposed heat flux, the flow can be choked 
within the tube. The minimally required 
temperature leading to ignition agrees quite well 
with the observed value. Simulations also suggest 
the need for measurements of the exit temperature 
and composition for accurate determination of the 
ignition condition. 
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